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ABSTRACT 



Aims. We present a catalog of sources of very high energy (E > 100 GeV) y-rays detected by Fermi telescope at Galactic latitudes 
\b\ > 10°. 

Metliods. We cross correlate the directions of individual photons with energies above 100 GeV detected by Fermi with the catalog 
of sources detected at lower energies. We find significant correlation between the arrival directions of the highest energy photons 
and positions of Fermi sources, with the possibility of chance coincidences at the level of 10"^^. We present a list of Fermi sources 
contributing to the correlation signal. A similar analysis is done for cross-correlation of the catalog of BL Lac objects with the highest 
energy photons detected by Fermi. 

Results. We produce a catalog of high Galactic latitude Fermi sources visible above 100 GeV. The catalog is split onto two parts. 
First part contains a list of 46 higher significance sources among which there can be 2 or 3 possible false detections. Second part of 
the catalog contains a list of 21 lower significance sources, among which 5 or 6 are possibly false detections. Finally we identify 7 
additional sources from the cross-correlation analysis with the BL Lac catalog. The reported sources of £ > 100 GeV y-rays span a 
broad range of redshifts, up to z ~ 1. Most of the sources are BL Lac type objects. Only 16 out of 74 objects in our list were previously 
reported as VHE y-ray sources. 

Key words. Catalogs - Gamma rays: galaxies - Galaxies: active - BL Lacertae objects: general 



1. Introduction 

Ground-based Cherenkov y-ray telescopes HESS, MAGIC and 
VERITAS have discovered a population of sources of Very 
High Energy (VHE) (E > 100 GeV) y-rays. Except for the 
sources discovered in the G alactic Plane survey by HESS 
jAharonian et aUl2005Ll2006l) . most of the sources were discov- 
ered in dedicated pointed observations. Surveys of large regions 
of the VHE y-ray sky with the existing Cherenkov telescopes 
are difficult because of the too narrow size of the field of view. 
Wide field of view ground-ba sed y-ray detectors MIL AGRO 
dAtkins et al.L |2004 and Tibet dAmenomori et all l2005h arrays 
have produced a systematic survey of the VHE y-ray sky. The 
energy threshold of the air shower arrays like MILAGRO and 
TIBET is rather high (in the multi-TeV band) so that only 
sources with spectra extending well above 1 TeV could be de- 
tected. 

All-sky monitoring of y-ray sources at the energies E ~ 
100 GeV became possible with the start of operation of Fermi 
telescope. Space-based Fermi telescope has a much smaller ef- 
fective collection area, compared to the ground based y-ray tele- 
scopes (~ 1 nr, compared to ~ 10'' m^ for the ground-based y- 
ray telescopes). At the same time, at the energies above 100 GeV 
Fermi signal is almost background-free (contrary to the ground 
based telescopes in which the signal has to be identified on top of 
strong backgroun d created by cosmic ra ys and optical/UV night 
sky background). iNeronov et al.l ( 1201 Ol) have searched for the 
point sources of £ > 100 GeV y-rays using Fermi data and found 



8 significant excesses with at least 3 photons within 0.1° circle 
corresponding to the 68% containment radius of point spread 
function (PSF) of the Large Area Telescope (LAT) on board of 
Fermi. Seven excesses were associated with known VHE y-ray 
sources, while the remaining one was identified with head-tail 
radio galaxy IC 310. Detection of IC 3 10 in the VHE b and was 
later confirmed by MAGIC telescope dMariottiL l2010l) . Due to 
the moderate collection area of Fermi, only the brightest VHE 
y-ray sources were detected individually in the Fermi VHE y- 
ray sky survey. All other known VHE y-ray sources at Galactic 
latitudes \b\ > 10° gave < 2 photons within the LAT PSF cir- 
cle and could not be found from the analysis of the data above 
100 GeV alone. 

A complementary method to identify the sources of £ > 
100 GeV photons detected by Fermi is to use a prior knowl- 
edge of source positions on the sky and to verify which of the 
already known sources could have produced the highest energy 
y-rays detected by the LAT telescope. In other words, sources 
of £■ > 100 GeV y-rays could be identified also via cross- 
correlation of arrival directions of the £ > 100 GeV y-rays with 
the source positions on the sky. Such an approach to the iden- 
tification of the sources was prev iously applied to the ana lysis 
of EGRET data above 10 GeV bv lDinpus & BretschI (1200 lb and 
iGorbunov etall (l2005h . 

In what follows we perform the cross-correlation analysis of 
the arrival directions of y-rays with energies above 100 GeV de- 
tected by Fermi at Gal actic latitudes \b\ > 10° with the 1-st year 
Fermi source catalog dAbdo et all 1201 Oah . Our analysis results 
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in a catalog of 46 high Galactic latitude sources of £ > 100 GeV 
y-rays. Although each individual source in the catalog is de- 
tected with significance around 3cr in this energy band, overall 
significance of detection of the entire source set is very high. 
Simple analysis shows that most of the sources contributing to 
the correlation signal are real VHE y-ray sources, only 2 or 3 
of them are expected to be false detections. Taking into account 
the fact that most of the sources in the list are BL Lac type ob- 
jects, we extend our cross-correlation analysis to the catalog of 
BL Lacs dVeron-Cettv & VeronL I20T0I) and find 7 more sources 
which correlate with the arrival directions of £ > 100 GeV y- 
rays and are not listed in the 1-st year Fermi catalog. Finally, we 
list for completeness sources from the 1-st year Fermi catalog 
for which the E > 100 GeV y-rays are found within the 95%, 
rather than 68% containment circle of LAT PSF. There are 21 
such sources, 5 or 6 of them are expected to be false detections 
due to the chance coincidence of the arrival direction of the VHE 
y-ray with the source position. 

The plan of the paper is like follows. In Section|2]we discuss 
data selection and data analysis methods. In Section[3]we present 
the results of the correlation analysis of the arrival directions of 
£ > 100 GeV y-rays with the sources of the 1-st year Fermi cat- 
alog. In Section|4]we give the list of sources contributing to the 
correlation signal. In Section [5]we perform the correlation a nal- 
ysis with the BL Lac catalog of lVeron-Cettv & Veronl ( l2010h and 
give the list of additional BL Lacs correlating with the highest 
energy Fermi photons, but not present in the 1-st year Fermi cat- 
alog. In Section|6]we comment upon individual sources. Finally, 
in Section|2]we discuss the results. 

2. Data selection and data analysis 

For our analysis we have used the LAT data which were col- 
lected in the period between August 4, 2008 and April 12, 2010. 
The data were filtered using gtselect tool, so that only y-ray 
events (evcls=3) with energies above 100 GeV were retained 
for the analysis. The resulting list of photons has 4427 events. 
A fraction of the £ > 100 GeV photons comes from the di- 
rections close to the Galactic Plane. Galactic Plane is a source 
of significant diffuse y-ray emission even at the energies higher 
than 100 GeV. Presence of strong diffuse emission complicates 
the analysis of the point source contribution. Taking this into ac- 
count we consider only photons coming from Galactic latitudes 
\b\ > 10°. The cut on the Galactic latitude leaves 2774 photons 
for the analysis. 

The set of photons considered in the anal ysis includes the list 
of photons studied bv lNeronov et alj (l2010h . In our analysis we 
ex clude the high-confid ence sources of VHE y-rays identified 
bv iNeronov et al. We also exclude 65 y-rays associated 

to these 8 sources from the Fermi £ > 100 GeV photon list. 
Otherwise, these sources would dominate the correlation signal. 
Nevertheless, we include the high-confidence sources in the final 
source list given in the Table [1] to obtain a complete catalog of 
sources of VHE y-rays found by Fermi. The final photon list 
used in the analysis includes 2709 photons with \b\ > 10°. 

Photons detected by FermijhPil are divided into two types 
- front- and ^acA:-converting. Photons that pair-convert in the 
top 12 layers of the tracker are classifie d as front-converting , in 
the other case they are back-converting ( lAbdo et al. L 1201 Obi) . At 
the same energy front-converting photons have somewhat better 
PSFQ- We take this fact into account and consider the sets of 

' http://www-glast.slac.stanford.edu/software/ 
IS/glas t_lat_performance.htm 



1028 front- and 1681 back- converted photons separately in our 
correlation analysis. 

If the sky region around the source is not too crowded, the 
point spread functions of individual Fermi sources do not over- 
lap at the energies above ~ 1 GeV. Besides, at high Galactic lati- 
tudes the strength and complexity of diffuse Galactic y-ray back- 
ground at the energies above 1 GeV is largely reduced (com- 
pared to the low-latitude / low-energy background). Taking into 
account these two facts we adopt a simplified procedure for the 
extraction of the source spectra. For each source and each energy 
we first calculate the number of source counts within a circle of 
the radius equal to the 95% containment circle of the LAT PSF 
and calculate the LAT exposure within this circle using gtexpo- 
sure tool. To estimate the diffuse sky background at the position 
of the source we calculate the number of background counts in 
a ring with the outer radius 3° and the inner radius equal to the 
95% containment radius of the LAT PSF, centered on the source. 
We verify if the region used for background estimate contains 
point y-ray sources (see e.g. the special case of the source RBS 
76 discussed in section |6]l. If a source from the 1-st year Fermi 
catalog, or a new source is present in the default background es- 
timate region, we either exclude a circle of the radius equal to the 
95% containment radius of the LAT PSF around this source from 
the background estimate region, or displace the background es- 
timate region altogether, so that it does not intersect with the 
95% containment radius of the LAT PSF around the "contami- 
nating" source. We have verified that the results obtained using 
such an intuitively simple spectral extraction procedure (similar 
to the one commonly used in X-ray astronomy) are consistent 
with those obtained via spectral extraction using likelihood anal- 
ysis in narrow energy bins and/or likelihood analysis in the wide 
energy range using a pre-defined broad-band spectral model. 



3. Correlation between arrival directions of 

E > 100 GeV photons and sources from the 1-st 
year Fenni catalog 
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Fig. 1. Dependence on the chance probability of correlation be- 
tween arrival directions of £ > 100 GeV photons and positions 
of Fermi sources as a function of the radius of the search cir- 
cle. Solid curve corresponds to all photons. Dashed and dotted 
curves correspond to respectively front and back converted pho- 
tons. 
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In order to find sources of £ > 100 GeV y-rays which in- 
dividually cannot be recognized in the Fermi data in this en- 
er gy band, we apply a me thod similar to the one di scussed 
bv lPingus & BretschI (1200 ih and iGorbunov et al.l (l2005h . Using 
this method, sources producing just a few photons could be 
identified. A fraction of identified sources could be false de- 
tections due to a chance coincidence between the source posi- 
tion and arrival direction of a background photon. The fraction 
of positive-to-false detections could be readily estimated. Real 
sources should largely outnumber the false source detections. 
This makes the list of sources contributing to the correlation a 
valuable "input catalog" for the observations with ground based 
y-ray telescopes which are more sensitive than Fermi in the VHE 
y-ray band, but, contrary to Fermi, do not have all-sky survey ca- 
pabilities. 

As a first choice, we took the 1-st yew: Fermi catalog as a n 
input catalog for the correlation analysis (lAbdo et al l l2010ah . 
The catalog contains 1451 sources, 1024 of which are objects 
with 1^1 >10 de grees. Removing 7 co nfirmed VHE y-ray sources 
from the list of iNeronov et al.l (1201 Ol) (IC 3 10 is not in the Fermi 
catalog), we obtain an input catalog of A^soui 



1017 sources. 




e [deg] 

Fig. 2. Number of events which come from Fermi sources as a 
function of the radius of the search circle (red data points). Error 
bars indicate fluctuations of background in determination of sig- 
nal from the sources. Average number of background photons is 
shown with green line. Vertical line at 0.1 degree is the radius 
of 68 % containment circle of Fermi PSF. Vertical line at 0.175 
degree shows the position of the dip in the probability shown in 
Fig- HI Horizontal line shows the asymptotic number of photons 
from all the sources. 



The probability for a y-ray to come within given distance 
from position of one source in the catalog is estimated as the 
ratio of the area of the circle of radius 9 and the part of the sky 
at 1^1 > 10°: 



Pi = 



47r(l -cos(80°)) 



^ 9.2 ■ 10"^ 



(1) 



The total probability for a photon to come within angle 6 from 
position of one of the catalog sources is then 



■pi = 9.3 ■ 10" 



1017 



O.P 



Probability that K or more photons from Ny - 2709 y-rays come 
within the angle 9 from any of sources in the catalog by chance 
is given by binomial probability 



k=K 



(Ny-k)\k\ 



(3) 



(2) 



In Fig. [T] we plot this probability as a function of angle 9. 
The probability that y-rays with energies higher than 100 GeV 
come from the directions close to the positions of Fermi sources 
by chance is P < lO^^***, i.e. the chance coincidence hypothesis 
is firmly ruled out. 

The function P{9) has two minima at around 0. 1 degree and 
0. 175 degrees. This is explained by the fact that precision of de- 
termination of arrival direction for the front and back converted 
y-rays is different. To verify this we split the whole photon list 
into two parts corresponding to the front and back converted y- 
rays and plot the contributions of those photons separately. From 
Fig.[T]one can see that minimum at 0. 1 degree is due to front pho- 
tons, while minimum at 0.175 degrees is due to both front and 
back photons. 

In Fig. |2]we plot the number of source photons with arrival 
directions within the angle 9 from the catalog of Fermi sources 
and compare it to the expected number of background events 
within the same distance from the sources. To estimate the num- 
ber of background events we used two alternative methods. First, 
we estimated the expected background near each source count- 
ing all photons within 10 degrees from it and estimating the frac- 
tion of photons which come within the angle 6. This method can 
be used only for 9 < 0.5°. For 6 > 0.5° circles around Fermi 
sources cover significant fraction of the sky and overlap. In the 
second method we estimated the background by generating 2709 
photons in the sky outside of Galactic plane \b\ > 10°, following 
Fermi exposure, estimated using gtexposure tool, and averaging 
over many Monte Carlo simulations. Both methods gave similar 
results for small angles. 

Vertical lines in Fig. |2] show angular positions of the two 
minima of probability seen in Fig.[T] The minima correspond to 
the places where the steep rise of the number of signal events 
changes to a plateau. Asymptotic of the signal counts at large 9 
gives an estimate of the total number of y-rays contributing to 
the correlation. One can see that only ~ 66 out of 2709 events 
come from the known Fermi sources. Ns, o.i - 40 of these 66 
events are displaced by less than 0.1° from the catalog source 
position. In this circle one expects only Nb.q.i - 2.5 back- 
ground events. No more than 6 background photons are expected 
at the 95% confidence level and no more than 8 background pho- 
tons are expected at the 99.5% confidence level. The number of 
events contributing to the correlation within 9 - 0.175° circle is 
Ns, 0.175 - 66, while the expected number of background events 
is Nb. 0.175 - 7.7 (12 at 95% confidence level, 15 at 99% confi- 
dence level). 

4. Catalog of extragalactic VHE y-ray sources 

The list of sources contributing to the correlation signal within 
= 0.1° is given in Table [1] There are A^source, o .i =46-7 - 
39 suc h sources, if the sources from the list of iNeronov et al.l 
(1201 Oh are excluded. On average, for Ny - 2709, only p ■ Ny ^ 
2.5 photons could happen to be at the distance 9 < 0.1° from a 
catalog source by chance. This means that 2-3 out of 39 new 
VHE y-ray sources listed in Table [T] are expected to be false 
detections (and no more than 6[8] sources are false detections at 
the 95% [99.5%] confidence level). 



4 A. Neronov, D.Semikoz and le.Vovk: High Galactic latitude Fermi sources of y-rays with energies above 100 GeV. 





IFGL 


RA 


Dec 


30-100 


30-100 


-'''68 
100-300 


A/ ■ 

^ 'nun 
L'*95j 


P 

100-300 


Type* 


Mnmp>9 
1 dlliCZ 




1 

i 


TOO^'^ S 1Q91 


8 '^01 ^ 

O. JV 1 O 


1 IfSSO 

- 1 7. JUjU 


U 


8 
y.o 


1 0- 


15 


1 h 


1 
i 


9 1 
Z.J 


. 1 0" 
• lU 


-b 


TIT V 


RPiQ 7A 
fx / U 


f\] 

U.U 1 




TOOSzL 9ASS 


1 T 7^17 
1 J. / J J / 


94 0900 
- ZJ-x .y LiyXj 



u 








1f 
il 




9 a 

Z.O 


. 1 0" 
• lU 


-3 










TOl 1 zL09^ 

Jul iU.U-H-UZj 


1 7 S9S0 

1 / . JZJU 


40 1887 

-H-U. JOO / 


1 

1 


1 9 
1 .z 


. 1 0" 
• lU 


_2 


1f 
il 




1 a 
1.0 


. 1 0" 
■ lU 


-3 








A 


10900 1 S990 

JuZUy. JT- JZZ7 


^9 1^09 
jz. J jyz 


S9 4007 

- jZ.H-yU / 


z 


8 1 
0. 1 


1 O" 


-5 


1f 
il 




9 
z.y 


1 O" 


-3 


TIT Vf 

L V 1 


R7Ti T0900 S990 

DZ^D JUZU" jzzy 






T091^ 9+9944 

J UZ i .ZTZZt-t- 


^99S 

J J . jZZ J 


22.7489 


2 


1 9 

1 .Z 


. 1 0" 
• lU 


-4 


lb 




9 9 
z.z 


. 1 0" 
■ lU 


-3 




C^ ASS T09 19+9944 

V-.J— irtOO JUZIZTZZt-t- 




f, 


T0999 ^^+4^09 

JUZZZ.UT^H-jUZ 


J J -UUO 1 


41 018S 

H-J.U JO J 


1 4 

IH- 


1 1 
J. 1 


1 0- 


35 


JL 1 Lf 


1 
i 


1 1 
J. 1 


1 0- 


12 


TIT V 

Di_, V 




444 


7 


T09^7 S '^f\(\'^ 


^0 17SS 

Dy .C 1 00 


lA OS 10 
- jU.Uj JU 


1 

J 


1 8 

1.0 


1 O" 


-7 


1 h 
i L) 




9 1 
Z. 1 


1 O" 


-3 


TIT V 


RRS 114 

fv . J JJH- 




Q 
O 


TO^O'^ S 940/1 

JUJUJ>. J-ZH-UU 


4S 8800 


94 1 1 94 

-ZH-. 1 IZH- 


8 




4 a 


1 o~ 
lU 


21 


1f 
il 


[il 


1 7 
1 . / 


1 O" 


-4 


TIT Vf 

DJ_i V 1 


PK'C 0101 941 
rivo ujUI ZH-j 


9(S 

u.zu 


Q 

y 


TO^I lfxC\Q 


48 OQH 


if\ 1 sa 1 

-ZD. 1 JO 1 


1 
1 


1 9 

l.Z 


• 10 


-2 






9 S 
Z.J 


1 O" 


-6 


TIT V 

V 


R V TOl 1 A 9 9A07 
IWV JUjlU.Z-ZDU/ 


441 
U.44j 


in 




40 OQOA 
^y.\jyy\j 


-D'+.U'+ 1 1 


n 
u 








i u 




1 9 
i.Z 


1 O" 


-3 








1 1 

1 1 


J UJZ J. 7- 


SI 4780 

J 1 .H- / OU 


1 A 8 1 74 

- lU.0 1 / H- 



u 








1 h 
i L) 




9 1 
Z. 1 


1 O" 


-3 


TIT V 


R R S 49 1 

fx IJ k J H-Zl 


90 

U.Z7 


1 ? 


TOl^S Oil 11 1 


S4 79 1 9 


11 9114 

1 J.Z J IH- 



u 








1 h 
i L) 




9 

z.y 


1 O" 


-3 








1 ^ 


T041 S-i-01 07 


A4 90/^0 


1 1 944 

1 . 1 ZH-H- 


n 
u 








1 f 
il 




9 9 
z.z 


1 O" 


-3 


TIT Vf 

D L V 1 


9F ndid-i-nn^7 

it/ utAt-ruus / 






JUH-ZO.U-J / JU 


A7 1 SfS7 

U / . 1 JU / 


17 041 9 
-J / .y-r 1 z 


1 

J 


S 9 
J.Z 


• 10 


-7 


1f 
il 




9 4 


. 1 0" 
• lU 


-3 


TIT 


PK'C 049 A 180 

l^JVo UH-ZU JOU 


1 1 1 

1.11 


1 ^ 


jujuj.y tuizi 


/ U.H-O J J 


f\\ 1S97 

Ul . J JZ / 


1 

1 


1 8 

1.0 


• 10 


-2 


1f 
il 




a 1 

O.J 


1 O" 


-7 








lU 


T0S07 0-i-/S7^S 


lf\ OQ^l 

/ \j.yyO 1 


U / .U JU / 


1 1 

1 1 


• 


10-2^ 


li 1 Lf 


i 


S A 
J.U 


1 O" 


-6 


TIT V 




141 

U. JH- 1 


17 
1 / 




SS OS 1 9 

oj.yj IZ 


SS SlOl 

-J J. J jUI 


1 
1 


1 A 

l.D 


1 r\' 


-2 


i U 




1 1 
J. i 


1 O" 


-3 


TIT V 


RyR T0S41 SS19 
0/. JUJ4j-JjjZ 




1 


TOf^SO 7-i-9S0^ 


1 09 fS81 7 

lUZ.UO 1 / 


Z J.UUZZ 


1 
1 


8 fi 
o.D 


. 1 r\- 


-3 


1f 
il 


n 1 
[il 


1 1 

1 . 1 


. 1 0" 
• lU 


-4 


TIT V 


1 FS 0A474-9S0 

lHjij UUH- / T^ZJU 


901 
u.zu J 


1 Q 


10710 /^-i-SOl 1 


1 07 /SSOS 

lU / -UJ" J 


so 1 asi 

jy. 1 J J 


1 

J 


1 .J 


. 1 r\- 


-6 


1f 
il 




1 A 
i.u 


1 O" 


-3 


TIT 




1 9S 

U. IZJ 




10791 0-4-7190 
ju/zi.yT^ / izu 


1 1 4704 

1 lU.H- / ^H- 


71 1448 

/ 1 . JH-H-O 


u 


a s 

O.J 


lU 


14 


1f 
il 




1 1 
J. 1 


1 O" 


-3 


TIT Vf 

L V 1 




1 

U.J 


9 1 


1074^^ 9a-74^S 
JU /H-J .Z-r / JO 


1 1 A ^9 1 S 
1 lU. JZ 1 J 


74 A4 1 S 


n 
u 








i u 




1 4 


1 O" 


-3 


TIT V 

J3J_i V 


1 FS 0717-1- 74A 
ICo U / J / /4U 


011s 

U. J 1 J 


99 


TOSOO 44-^4SS 

JUOUy.H-T- JH-J J 


1 99 ^S'^Q 

1 ZZ. J J J7 


14 Oil 1 

jH-.y J J 1 



u 








1f 

il 




9 4 


. 1 0" 
■ lU 


-3 


TIT Vf 

DJ_i V 1 


R9 080A-1-1S 

DZ UOUUtjJ 


089 S 

U.UOZ J 


9"^ 


TOSOO S-4-S9 1 
juouy.jTjziy 


1 99 ^8S^S 

1 ZZ. JO JU 


S9 1178 

jZ. J 1/0 


1 

J 





• 10 


-7 


1 f 1 h 

li 1 L) 




A A 


1 O" 


-6 


TIT Vf 

L V 1 




1 18 

U. 1 JO 


94 


TOSl 4 1^11 

JUO i U.H— 1 J 1 i 


1 94 1 1 9 1 

1 ZH-. 1 IZ 1 


1 1 1 oil 
- 1 J. 1 7 J J 


1 

J 




• 10 


-8 


1f 
il 




1 /I 


1 ri" 


-3 


TIT Vf 

DJ_i V 1 


R7R T081A 1111 

DZjD JUOlU ljll 




9'^ 


TOOO'^ '^-1- 1 


1 JU. jy ID 


1 1 0481 


1 
1 


1 1 


• 10 


-2 


i u 




1 1 
J.J 


1 ri" 


-3 




r'R ATFS TOQOSj-IISa 
\^Krtl J_<ij JU^UJt^IjJo 




9A 
ZD 


T001 7-1-90^ 1 
j\}y lJ.I tZ7 j 1 


1 18 041 9 
1 jo.y'+lZ 


90 SI 1 


1 

J 


S 1 
J.J 


• 10 


-7 


1 f 

il 




9 1 
Z. 1 


1 O" 


-3 


TIT Vf 
L V 1 


R9 001 9-u90 
DZ U^IZt-Z^ 




97 

Z / 


JUyj j.U-Uojo 


1 48 lf\l'\ 
I'+o.Zu / J 


8 A4A1 


7 

z 


u. / 


• 10 


-5 


1 f 

il 




A 
Z.4 


1 r\- 


-3 


TIT Vf 
D L V 1 


RVS T00S10 0840 
rwvo JU7JJU-U04U 




98 
Zo 


T00S7 7_i-S'>91 
JUyj / . / -rj jZj 


1 40 411 S 
I'+y.f jl J 


SS 1099 

jj. jyzz 


A 
H 


n 1 
/.J 


• 10 


-9 


1 f 

il 




1 1 
J.l 


1 ri" 


-3 




AC -u S S 17 
4\w i- J J . 1 / 


aoss 

U.oVJ J 


90 

zy 


T1 01 '^ 1 -j_4097 


1 SI 788S 
1 J J. / OOJ 


4Q 4S9A 


J 


1 /I 


10- 


11 


1 f 


1 
1 


1 A 
J. 4 


1 r\- 


-5 


T!T \/ 
DL, V 




919 
U.ZiZ 


jU 


T1 1 04 4-1.^81 9 
J 1 IUt-.t-t^ JO 1 Z 


1 fxfx 1 947 
1 UU. IZH- / 


18 9 1 08 
JO.ZIUO 


J 1 









i Ji / L> 






U 




TIT V 


iVARIl *tZA 


01 
U.Uj 


D 1 


11117 1 4-901 

Jill / . 1 T^ZU 1 J 


1 f\Q 900"^ 
1 \jy .ZJyy J 


90 9999 

Z.\j.AZ,Z.Z. 


1 

J 


1 7 
J. / 


• 10 


-7 


1f 
il 




9 S 

Z.J 


1 O" 


-3 


TIT Vf 

Di_, V 1 


PR ATFS Tl 1 1 74-901 4 

\^J\./\l CjO Jill/ T^ZUIH- 


1 

U. 1 


^9 
JZ 


T1 9 1 7 7-1-^007 
J 1 Z 1 / . / -r jUU / 


1 84 AAf>rK 
1 OH-.*+'+D J 


10 1 90A 

jU. 1 zuu 


7 

z 


1 A 


• 10 


-4 


1 f 
il 




9 S 
Z.J 


1 O" 


-3 


TIT V 


R9 1 9 1 S-i-10 
DZ IZlj-rjU 


1 1 
U. 1 J 




T1 994 7-J-9 1 9 1 
J 1 ZZH". / -rZ IZl 


1 8A 1 087 

1 OU. lyo / 


91 1A11 

Z 1 . JU J J 


u 








1 f 
il 




1 7 


1 ri" 


-3 




AC -1-9 1 IS 
4\w. i-Zl.jJ 


419 
U.4jZ 




Tl '^^7 7 1 9SS 

J i-JJ 1 . /-IZJJ 


904 AXXfx 

ZUH-.H-J jU 


1 9 0980 
- iz.yzo? 



u 








1f 
il 




1 1 
1 .J 


1 O" 


-3 




PK'C 1 lis 1 97 

l^JVO IJJj-lZ/ 


S10 

U.J J" 




Tl 4^7 0-i-S^^40 
J If J / .U+jD'+U 


9 1 9A1 7 
ziy.zui / 


SA A790 
jD.O / ZU 


1 

J 


8 1 
O.J 


• 10 


-7 


1 f 
il 




1 
J.U 


1 O" 


-3 


TIT 


RyR TI4IA4-SAIO 

DZjD J i4JOT-JUJ7 


1 s 

U. 1 J 


jD 


Tl 'iS'^ 7-1- 1111 

J 1 J J J . / T 1 1 1 1 


918 0181 

ZJ0.7JO J 


11 1 09 1 
1 1 . lyz 1 


9S 

Zj 









ZiZU 


1 
i 


9 7 
Z. / 


1 n~ 

lU 


12 


TIT V 
Di-. V 


PP 1 ^^'^4-1 1 ^ 
r vj A33J-rAAj 




J / 


^^f\'^'^ 04-^04s 

J luj j.vt^ jyT-j 


9S1 4807 

ZJ J.H-O" / 


10 7S97 
Jy. 1 jZ / 


1 4 

IH- 


7 9 
/ .Z 


10- 


37 


Of 

71 


1 
i 


A 


1 0- 

lU 


29 


TIT V 


iVAti.ll Z)UA 


OllfS 
U.UJ ju 


Do 


T1799 S-i-1019 

Jl / ZZ.jTlUlZ 


7^0 /i404 

ZUU.UH-yH- 


10 9 1 44 

lU.ZlH-H- 



u 








1f 

il 




9 9 
Z.Z 


1 O" 


-3 




TRATFS T17994-1011 

X^JVrtl CjO J L / ZZt^ lU 1 J 




"^Q 
jy 


Tl 744 9-1- 1 0^4 
J 1 /H-'+.Z-r l^jH- 


lf\f\ 071 
ZUU.U / lU 


1 S70A 
1 7. J /UU 


1 
1 


1.9 


• 10 


_2 


i u 




2.2 


. 10" 


-3 


TIT V 


LlLo l/41i^l7U 


081 
U.Uo J 




Tl 894 f\A- 101 

J 1 OZH-.Ut^ lU 1 J 


97A 1 71Q 

Z / U. 1 / J7 


1 9 1 80 

LKJ.Z, 1 OU 



u 








1f 

il 




1.7 


■ 10" 


-3 








41 


J2000.0+6508 


300.0215 


65.1334 


5 


5.4 


10- 


11 


If 


ni 

L ^ 1 


1.4 


• 10" 


-4 


BL Vf 


lES 1959+650 


0.048 


42 


J2004.8+7004 


301.2074 


70.0710 











lb 




3.8 


• 10" 


-3 








43 


J2009.5-4849 


302.3761 


-48.8295 


7 


4.6 


10- 


16 


4f 


[1] 


2.4 


10" 


15 


BL V 


PKS 2005-489 


0.071 


44 


J2158.8-3013 


329.7117 


-30.2176 


22 









5f3b 


1[1] 


1.2 


10- 


24 


BL V 


PKS 2155-304 


0.117 


45 


J2314.1 + 1444 


348.5487 


14.7463 


1 


1.6 


• 10 


-2 


If 




2.4 


■ 10" 


-3 


AGN 


BZU J2313+1444 


1.319 


46 


J2322.6+3435 


350.6540 


34.5839 


1 


1.9 


• 10 


-2 


If 




3.0 


• 10" 


-3 


BL 


CRATES J2322+3436 


0.098 



Table 1. List of Fermi sources for which a £ > 100 GeV photon is found within the 68% containment circle of radius 0.1°. 1-st 
column: Fermi source ID; 2-nd column: Right Accention; 3-rd column: Declination; 4-th column: number of photons in the 68% 
containment circle in 30-100 GeV band; 5-th column: chance coincidence probability to find given number of photons in the 68% 
containment circle in 30-100 GeV band; 6-th column: number of photons in the 68% containment circle in 100-300 GeV band; "f 
marks front photons, "b" marks back photons; 7-th column: number of photons in the < 0.175° circle around the source or in the 
95% containment circle of the radius 0.3° (square brackets); 8-th column: chance probability to find given number of background 
photons within 0.1° and 0.175° [0.3°] circles around the source position; 9-th column: alternative name of the source; 10-th column: 
source redshift. Bold marks the sources which have been previously detected in the VHE band by ground-based y-ray telescopes. 
' BL - BL Lac type object; QSO - quasi-stellar object, AGN - active galactic nucleus, FSRQ - flat spectrum radio quasar; FR I - Fanaroff-Riley 
type I radio galaxy; V - also found from the correlation analysis with the BL Lac catalog; Vf - front photon contributes to the correlation signal 
with the BL Lac catalog. 



It is not possible to know in advance which of the 39 sources 
from Table [1] are false detections. In principle, the 8 sources 
which have more than one photon within the 95% containment 
circle of the LAT PSF are not expected to be false detections, be- 
cause the probability that 2 background photons come close to 



one and the same source is much smaller than the same proba- 
bility for just one background photon. Additional information 
which might help to single out the false detections could be 
obtained from comparison of the estimate of the source fluxes 
expected from the extrapolation of the measured spectral char- 
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Fig. 3. Comparison of the distribution of fluxes and photon in- 
dices of Fermi sources of E > 100 GeV y-rays (triangles) with 
those of all sources from the Fermi AGN catalog. Blue triangles 
mark known TeV y-ray sources. 



acteristics of the sources in the 1-100 GeV energy band with 
the simple estimates of the flux which produces 1 photon at 
E > 100 GeV within 1.5 year exposure with Fermi. In princi- 
ple, a false detection might be spotted if the source contributing 
to the correlation has low flux and/or soft spectrum. This is il- 
lustrated by Fig.[3]in which the distribution of fluxe s and photon 
indices of the AGN from the Fermi AGN catalog dAbdo et al.L 
1201 Od) is compared to the distribution of fluxes and photon in- 
dexes of the sources from Table [1] Dashed lines show the com- 
binations of the spectral parameters which are expected to give 
fixed number of photons in the 100-300 GeV band, assuming no 
high-energy cut-off in the spectrum. 

One can see that the "candidates" for the false detections 
are the sources with softest spectra, which are 4C -H21.35, PKS 
1335-127 and CRATES J1722+1013. This is also clear fromFig. 
|4]in which the measurements of the source fluxes in 0.3-1 GeV, 
1-3 GeV, 3-10 GeV and 10-100 GeV from the 1-st year Fermi 
catalog are compared with the estimates of the source fluxes in 
the 100-300 GeV band based on the number of £ > 100 GeV 
detected photons from each source. 

Although the significance of detection of individual sources 
of 1(2) photon in the 100-300 GeV band in Table [1] is around 
3cr(4cr), some of the sources are detected with significance 
higher than 5cr in the adjacent 30-100 GeV energy band. The 5- 
th and 6-th columns of Table[T]show the numbers of photons and 
the chance coincidence probabilities derived from comparison 
of the expected number of background photons within < 0. 1 2° 
circles (corresponding to the 68% containment circle at 30 GeV) 
around the sources with the number of source photons in the 30- 
100 GeV band. For the sources which are detected with more 
than 5cr significance in the 30-100 GeV band we show more 
detailed spectra above 1 GeV, calculated using all publicly avail- 
able data of Fermi collected from August 2008 to April 2010 
(rather than 1 1 month data used for the analysis of the 1-st year 
Fermi catalog). One can see that for most of the sources the es- 
timates of the source flux in the 100-300 GeV energy bin agree 
well with the extrapolation of the powerlaw spectrum from lower 
energies. 

Maximum of the coiTelation signal (minimum of the chance 
probability) shown in Fig.|2]is achieved at the angle - 0.175°, 
corresponding to the 68% containment radius for the back type 



events. Improvement of the correlation with the increase of 
from 0.1° to 0.175° means that there are more real events which 
came within the rings 0.1 < < 0.175° around their sources, 
than there are background events in these rings. 

Table |2] lists the sources contributing to the correlation sig- 
nal within 6 < 0.175° but not within 6 < 0.1°. The number of 
background events in the ring 0.1° < 9 < 0.175° is estimated to 
be 5.4. This implies that 5 or 6 out of the 21 additional sources 
listed in Table|2]might be false detections. 

It is important to note that only 16 out of 67 sources listed 
in Tables [T] and |2] are akeady known sources of VHE y-rays 
detected by the ground-based y-ray telescopes. Majority of the 
sources from Tables [1] and [2] are new real detections in the VHE 
band. Taking into account the fact that only 8 out of the 67 
sources listed in Tables [1] and |2] are possible false detections, 
Fermi all sky survey at the energy Ey > 100 GeV reveals 67-8- 
16=43 new extragalactic VHE y-ray sources. This doubles the 
number of already known extragalactic VHE y-ray sources 0. 

Another important point is that the E > 100 GeV y-ray 
sources listed in Tables[T]and|2]are distributed over a broad range 
of redshifts. The "record" redshift of the VHE y-ray source so 
far was z - 0.5 for the possible de tection of 3C 279 d uring a 
short several hour flare by MAGIC dAlbertet al.Ll2008h . In our 
Table [T] several sources have redshifts larger than 0.5. 



5. Correlation between arrival directions of 

E > 100 GeV photons and sources from the 13th 
Veron catalog. 




1e-40 



Fig. 5. Dependence on the chance probability of correlation be- 
tween arrival dir ections of £ > 100 GeV pho tons and positions 
of BL Lacs from lVeron-Cettv & VeronI (1201 Ol) catalog as a func- 
tion of the radius of the search circle. Notations are the same as 
in Fig.[T] 

Most of the £ > 100 GeV sources listed in Tables [T] and |2] 
are BL Lac objects. These sources are known to be c haracterized 
by rela tively hard spectra in the GeV energy band (lAbdo et al.L 
l2010d) . It is, therefore, possible that some of the VHE y-ray 
emitting BL Lacs escape detection in the GeV band, while being 
bright VHE y-ray sources. 



See |http://tevcat.uchicago.edii71 
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1FGL J0033.5-1921, spec. Index = 1.89 | 1 FGL J0054.9-2455, spec, index = 1.95 



1FGL J01 10.0-402J, spec, index = 1.34 



IFGL J0209.3-5229, spec. Index = 2.01 



1FGL J0213.2 + 2244, spec, index = 1 



1FGL J0222. 6+4302, spec, index = 1.93 



IFGL J0237. 5-3603, spec, index = 1.21 



IFGL J0303.5-2406, spec, index = 1.98 



IFGL J0315.9-2609, spec, index = 1 



IFGL J0316.3-6438, spec, index = 2.03 



IFGL J0325.9-1649, spec, index = 1,88 



IFGL J033a,8+1313, spec, index = 1.78 



IFGL J0416.S + 0107, spec, index = 1.3 



IFGL J0428.6-3756, spec. Index = 2.13 




IFGL J0505.9 + 6121, spec, index = 1.73 



i i 



IFGL J0507.9 + 673a, spec. Index = 1.75 



IFGL J0543.S-5531, spec, index = 2.05 | IFGL J0650.7 + 2503, spec, index = 2.04 



IFGL J07 10.6+591 1 , spec, index = 1.28 



IFGL J0721,9 + 7120, spec. Index = 2.15 



IFGL J0745.2 + 7438, spec, index = 1,70 | IFGL J0809,4 + 3455, spec. Index = 1.77 



IFGL J0809.5+5219, spec, index = 2.10 



IFGL J0816,4-131 1 , spec. Index = 1.83 



IFGL J0905.5+1356, spec, index = 1.95 



IFGL J0915.7 + 2931, spec. Index = 1.9' 



IFGL J0953.0-0838. spec. Index = 1.9C 



IFGL J0957.7 + 5523, spec. Index = 2.05 



IFGL J1015.1+4927, spec, index = 1,93 



IFGL J1104,4+3ai2, spec, i 



IFGL J1117.1+2013. spec, index = 1,71 



IFGL J1217.7 + 3007, spec. Index = 1 



IFGL J1224. 7+2121, spec, index = 2,55 



IFGL J 1337.7- 1255, spec, index = 2.50 



IFGL J1437.0+5640, spec, index = 1.45 



1 FGL J 1 555.7+ 1111, spec. Index = 1 .66 



IFGL J1653,9 + 3945, spec, index ^ 1,65 



IFGL J1722, 5+1012, spec, index ^ 2,33 



IFGL J 1744.2+ 1934, spec, index ^ l.e 



IFGL J1824, 6+1013. spec. Index ^ 2.06 



i i i i 



IFGL J2000,0 + 6506, spec, index = 2,10 



IFGL J2004,6 + 7004, spec, index = 1.94 



IFGL JZ009.5-4849, spec, index = 1.90 



IFGL JZ158-8-3013. spec. Index = 1.91 



IFGL J2314, 1 + 1444, spec. Index = 2.12 



IFGL J2322, 6 + 3435, spec. Index = 1.82 



10^ 10^ 10^ 



10^ 10^ 10^ 



Energy, MeV 

Fig. 4. Comparison of the spectra of the sources from Table [1] in the 0,3-100 GeV band with the estimated source fluxes in the 
100-300 GeV band (green triangles). Blue circles (blue arrows) are the flux measurements (upper limits) from the 1-st year Fermi 
catalog. Green boxes are re-calculations of the source spectra using 1,5 yr exposure data set. 
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Table 2. List of Fermi sources for which a E > 100 GeV photon is found within the radius - 0.175°, but not within 6 — 0.1°. 
Notations are the same as in Table 1 . 
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Fig. 6. Number of events which come from BL Lac sources as a 
function of the radius of the search circle. Notations are the same 
as in Fig. 121 



To search for the VHE y-ray emitting BL Lacs with very hard 
spectra we have performed the correlation analysis similar to the 
one reported above for the 1-st year Fermi catalog also for the 
13th Veron catalog of BL Lacs dVeron-Cettv & VeronlQoiOl) . In 
order to select a co mplete sample o f BL L acs, we have adopted 
the prescription of iGorbunov et al.) ( l2005l) and made a cut on 
magnitude of objects Y < 18.0, which leaves 316 BL Lacs. This 
set of objects includes 6 out of 7 known V HE y-ray sources iden - 
tified in the E > 100 GeV Fermi data bv lNeronov et all (l2010h . 
Similarly to the analysis of the 1-st year Fermi catalog, we re- 
moved those objects and associated photons from analysis. We 
have also removed the 3 photons associated with IC 310. This 
leaves A^source = 310 BL Lac objects and 2709 photons. 

The probability for a single photon to come by chance near 
one of the selected BL Lacs is p = Nsoum • Pi - 2.7 ■ 10-"* ■ 



(6/O.rr. The probability that K or more photons fall within 
an angle from one of the sources is calculated using Eq. [3] 
Fig. |5] shows this chance probability as function of 0. Similarly 
to the case of the analysis of 1-st year Fermi source catalog, the 
probability has a minimum at the angular resolution of LAT at 
the energies £ > 100 GeV. 43 BL Lac out of 310 correlate with 
Fermi data within the angle < 0.175° corresponding to the 
minimum of P{0) from Fig.[T] 



In Fig. |6] we compare the number of source photons to the 
number of background photons as functions of the angular dis- 
tance to the BL Lacs. Asymptotical value of the signal events 

= 44 is reached at the angular distance ~ 0.2°. At the an- 
gle 0.1° 26 photons come from 25 sources, while only 0.77 of 
them are in average due to background. At 95 % C.L. there are 3 
background photons and 4 at 99% confidence level. At the angu- 
lar distance 0.175° 41 photons come from 35 sources, and 2.4 of 
them due to background in average. There are 6[7] background 
photons at the 95%[99%] confidence level. 7 out of those 35 
sources are not present in Tables[T]and[2] These seven sources are 
listed in the Table[3] Other sources which contribute to the corre- 
lation signal with both 1-st year Fermi and lVeron-Cettv & VeronI 
(l20Tol) catalogs and marked as "V" in Tables [T| and |2l 



We have also performed special analysis with only front 
photons and BL Lac catalog. In this case the correlation sig- 
nal is maximal at ^ 0.07° (see Fig.|5]). At this angle there are 
16 forward photons, which correlate with 16 BL Lac sources. 
Only 0.14 background photons are expected in this angle, which 
means that all the 16 sources are most probably real. At the 99% 
confidence level 2 of the 16 sources might be false detections. 12 
of these 16 sources are marked as "Vf in Tables[T](l 1 sources) 
and|2](l source) and the remaining 4 are those which have front 
photons marked "f in Table[3] 



8 



A. Neronov, D.Semikoz and le.Vovk: High Galactic latitude Fermi sources of y-rays with energies above 100 GeV. 





Name 


D A 

KA 


Dec 


A! 

-"'0.175 


V 
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lb 


0.071 
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B3 0651+428 
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lb 
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3 


RXS J09130-2103 
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Iflb 


0.198 


4 


RXS J10162+4108 
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lb 


0.281 


5 


SDSSJl 14023.48 


175.098 


15.469 


If 


0.244 




+ 152809.7 
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MS 12218+2452 


186.101 


24.607 


If 


0.218 


7 


MS 13121-4221 


198.764 


-42.614 


If 


0.105 



catalog for which a £ > 100 GeV photon is found within 
the circle of radius 0.175° and which are not the 1-st year Fermi 
catalog sources. 
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Fig. 7. Broad band spectrum of RBS 76. Grey data points are his- 
torical data from NED. Blue thin data points show the spectrum 
of an additional source FERMI J0535.7-1849 visible in Fig. [8] 

6. Comments on individual sources 

6.1. High redshiU sources 

Several sources in Table [T] have redshifts higher than the red- 
shift of 3C 279, z - 0.5, which i s the furtherest known 
VHE y-ray source (lAlbert et all l2008h . Study of the high red- 
shift sources in the VHE band is interesting for the mea- 
surements of the density and cosmologic al evolution of the 
Extragalactic Background Light (EBL ) ("Gould & Schredei , 
1967HKneiske et al.i l2004t IStecker et al.L l200(k iMazin & Raud 



2007 : iFranceschini et al.L l2008h and of the cosmological 



magn etic fields (iNeronov & Semikozl l2009t iNeronov & VovkL 
I2OIOI) . In this subsection we provide some details about the high 
redshift sources from Table [1] 

RBS 76 atz = 0.61 has 2 photons within the distance 0.175° from 
the catalog source position. It is detected with significance more 
than 4cr in the 100-300 GeV band. The source is also detected 
with significance more than 7cr in the 30-100 GeV energy band. 

Broad band spectrum of the source is shown in Fig. |2l 
The historical multiwavelength data are taken from the NASA 
Extragalactic database (NED) 0. Comparing the radio. X-ray 
and y-ray data one could conclude that the broad-band spectrum 
of the source fits into the general scheme of synchrotron-self- 
Compton or synchrotron-external Compton models with radio- 
to-X-ray emission being produced via synchrotron mechanism 
and y-ray flux being generated via inverse Compton scattering. 

In the course of investigation of the sky region around 
RBS 76 in the 1-300 GeV energy band we have discov- 
ered an additional source which is not reported in the 1- 
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Fig. 8. Fermi count maps of the sky region around RBS 76 in 
1-10 (left) and 10-100 GeV (right) energy bands. Circle shows 
the new unidentified Fermi source. Coordinate grid {RA, DEC) 
spacing is 1°. 
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Fig. 9. Broad band spectrum of PKS 0426-380. Grey data points 
show historical data taken from NED. 



St year Fermi catalog. The new source FERMI J0535.7- 
1849 has coordinates RA=5.56 + 0.13 DEC=-18.86 ± 
0.12. We have performed standard Fermi likelihood analy- 
sis (http : //fermi . gsf c . nasa . gov/ssc/data/analysi s/) 
and found that the Test Statistics (TS) value which character- 
izes the likelihood of source detection is TS= 136.5, which cor- 
responds to the sou rce detection at approximately llcr level 
dMattox et al.L Il996l) . There are no obvious low energy coun- 
terparts of the source. The source spectrum, shown in Fig.|7]is 
peaked in the 30-100 GeV band. The spectral characteristics of 
the source and the absence of low energy counterparts are con- 
sistent with the possibil ity that the source is a "blazar remnant" 
of the type discussed bv INeronov et al.l ( 1201 Obi) . 

PKS 0426-380 at z - 1 . 1 1 is a bright source in the GeV energy 
band, with the flux reaching 10"'" erg/cm^s. Its y-ray spectrum 
is relatively soft and is consistent either with a powerlaw with 
photon index F > 2.5 or with a cut-ofF powerlaw type spectrum. 
The y-ray measurements are compared to the multiwavelength 
data on the source in Fig.|9] Similarly to the case of RBS 76, the 
overall spectrum of the source could be interpreted in terms of 
the synchrotron self-Compton or synchrotron external-Compton 
type models with radio-to-optical emission being produced via 
the synchrotron mechanism and X-ray-to-y-ray emission being 
produced via inverse Compton scattering. Extremely high red- 
shift of the source and possible presence of the high-energy cut- 
off in the spectrum might make detection of the source with 
ground-based y-ray telescopes at the energies significantly larger 
than 100 GeV difficult. 



A. Neronov, D.Semikoz and le.Vovk: High Galactic latitude Fermi sources of y-rays with energies above 100 GeV. 



-10 



-12 



_ AC ■ 55.17 {Z-0.B955) 



EGRET 

^ Fermi 



-12 



_ PKS 0301 243 (z-0.26) 




Fcriiii _ 




■ 




■ 






* 









-5 



10 



-5 



10 



log K [eV] 



log K [eV] 



Fig. 10. Broad band spectrum of 4C +55.17. Grey data points Fig. 12. Broad band spectrum of PKS 0301-243. Grey color 
show historical data from NE D. Blue data point is the measure- shows historical data from NED. 
ment of y-ray flux by EGRET ([Thompson etal.Lll995h . 

n I I I I I I I I I I I I r 
BZU J2313 I 1444 (z-1.319) 



6 

-i; -12 



Fermi 



-14 



- 1 ' ' 

_ RX J0316.2 


1 1 

2607 (z-0.443) 


1 1 


Fcrnii _ 




* 
« 




+ 




• 

1 


1 





-5 




log 1 


5 

[eV] 


10 



^''^^ Fig. 13. Broad band spectrum of RX J03 16.2-2607. Grey color 

Fig. 11. Broad band spectrum of BZU J2313+1444. Grey data shows historical data from NED. 
points show data from NED. 



4C +55.17 SLlz — 0.8955 also has a relatively soft spectrum in the 
1-300 GeV band. The y-ray measurements are compared to the 
broad band data in Fig. [10] The source was previously detected 
in the 0.1-1 GeV band by EGRET. Comparison of the historical 
EGRET measurement with the Fermi measurement of the source 
flux shows that the source is variable at 10 yr time scale. 

BZU J2313+1444 is formally the highest redshift source in Table 
[1] with z - 1.319. The source is relatively weak. Contrary to the 
other high-redshift sources, the significance of source detection 
in the 30-100 GeV band is less than 5cr. From Fig.[TT]one can 
see that the source spectrum does not show a signature of high 
energy cut-off. This implies that the source might be an ideal 
candidate for the study of EBL evolution with redshift. 

6.2. Sources detected at > 4cr level atE> 100 GeV 

Ap art from the 8 h i gh-con fidence sources at £ > 100 GeV found 
by iNeronov et alj ( 1201 Ol) and the bright high-redshift sources 
discussed above, several other sources in tables [1] and |2] have 
detection significance higher than 4cr above 100 GeV. Below we 
list these sources and discuss some details of their broad band 
spectral energy distribution. 

PKS 0301-243 at z = 0.26 has one back and one front- 
converted photon associated with the source above 100 GeV. The 
back photon is within the 95% containment circle of the radius 
6 = 0.3°. Taken together the two photons result in a ^ 3.8cr 
detection of the source above 100 GeV. At the same time, the 



source is also detected with high significance in the 30-100 GeV 
band. The broad band spectrum of PKS 0301-243, shown in Fig. 
[T2l is consistent with the possibility that radio-to-X-ray emis- 
sion are produced via synchrotron mechanism and 1-300 GeV 
emission via the inverse Compton mechanism. 

RXJ03 16.2-2607 is a relatively high redshift (z = 0.443) source 
which has 2 backy-rays within the distance 9 = 0.1°. The chance 
coincidence probability of finding two photons within this dis- 
tance from the source is 2.5 x 10"* which implies detection sig- 
nificance 4.7cr above 100 GeV. From Fig. [T3]one can see that 
the source has hard y-ray spectrum without a signature of high- 
energy cut-off. The y-ray energy flux from the source is com- 
parable to the optical flux. Similarly to BZU J2313+1444, this 
source also might be an ideal candidate for the study of EBL and 
of cosmological magnetic fields. 

IFGL J0505.9+6121 has 3 photons next to it within 0.22°. 
Formally, the chance probability of finding three photons inside 
this distance is P = 8 ■ 10"*^, however, since detection angle was 
not defined a priori, it has to be penalized for this and the final 
probability is just 5cr (see Table [TJ, which still implies that the 
source is detected at E > 100 GeV individually, and not only 
within the population of Fermi sources from Table [T] It is also 
detected with significance higher than 5cr in the 30-100 GeV 
band. The spectrum of the source shown in Fig. [14] is consis- 
tent with a hard powerlaw (photon index harder than 2) up to the 
highest energy. 

The object is located near the Galactic Plane atb - +12°, so 
that it can be of Galactic origin. High level of Galactic diffuse 
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Fig. 14. Fermi spectrum of IFGL J0505.9+6121. 
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Fig. 15. Broad band spectrum of lES 0647+250. Grey color 
shows historical data from NED. 



y-ray background around the location of the source prevents its 
detection at lower energies E < I GeV. 

Detailed multiwavelength observations are needed to con- 
strain its nature. A possible candidate for the source identifi- 
cation is the radio- an d X-ray loud AGN RX J0505.9+6113 
(iBrinkman et al.L Il997h situated at the distance 7.6' from the 
catalog source position and at the distance 0.06° from the front 
£■ > 100 GeV photon. 

IBS 0647+250 at z = 0.201 was alre ady mentioned as a "can- 
didate " TeV blazar in the analysis of ICostamante & GhiseUinj 
(l200l . based on its broad band spectra properties. Fig. [15] 
demonstrates that the source has hard y-ray spectrum, with no 
signature of cut-off up to 300 GeV. This assures that it should be 
readily detectable with the ground based y-ray telescopes. 

6.3. Sources from the BL Lac catalog 

iJX709750-2705has 2 y-rays with energies £■ > 100 GeV within 
0.1° distance. One of the y-rays is a front type y-ray, for which 
the 68% containment circle has smaller radius. The front y-ray is 
at the distance 0.01° from the catalog source position and a back 
photon at 0.03°. The field of the radius 10° around the source 
contains 20 photons at the energies above 100 GeV. The overall 
chance probability of finding a front and back photons within the 
68% containment circles of the front and back photons is ^ 10"^. 

In spite of the fact that the source is not reported in the 1-st 
year Fermi catalog, it is detected with Fermi also at the ener- 
gies below 100 GeV. Standard likelihood analysis in the energy 
range 1-300 GeV gives the TS value 80 which corresponds to 
the significance of source detection ^ 9cr. 
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Fig. 16. Broad band spectrum of RX J0913.0-2103. Grey data 
points show historical data from NED. 
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Fig. 17. TS map of the region 2° x 2° in the energy band above 
1 GeV around the position of RX J0913.0-2103. 



The map of TS values in the region 2° x 2° around the source 
is shown in Fig. [17] The TS values were calculated for the model 
of source distribution in the region 14° around RX J09 130-2 103 
in which all the the 1-st year Fermi source catalog sources except 
for IFGL J0908. 7-21 19 were included and one additional source 
was allowed to have variable position. One can see that the TS 
value at the position of RX J09 1 30-2 1 03 is much larger than that 
at the position of IFGL J0908.7-21 19 which is only marginally 
detected above 1 GeV. 

The broad band spectrum of the source is shown in Fig. [16] 
One can see that the y-ray spectrum of the source is hard, which 
explains its non-detection in the first 1 1 months of operation of 
Fermi. 

MS 12218+2452 has a front photon at the distance 0.04 from 
the catalog source position. The source is also not reported in 
the 1 -St year Fermi catalog, but is detected in the 1 .5 yr exposure 
with TS value 80, which corresponds to the ^ 9cr significance of 
the source detection. 

The map of TS values in the region of the size 1° x 1° at the 
energies E > I GeV around the source position is shown in Fig. 
[T9l To construct this map we included all the sources from the 
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Fig. 18. Broad band spectrum of MS 12218+2452. Grey data 
points show historical data from NED. 
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Fig. 19. TS map of the region 1° x 1° in the energy band above 
1 GeV around the position of MS 12218+2452. 



1-st year Fermi catalog situated at the distance < 14° from the 
source position into the model of source distribution on the sky. 
One can see that the source can be clearly in the TS map. 

Fermi data are compared to the broad band archival data in 
Fig. [18] y-ray flux from the source is comparable to the opti- 
cal flux. The y-ray spectrum does not show signatures of high- 
energy cut-off. 

7. Discussion and conclusions 

Study of cross correlation of the arrival directions of the highest 
energy photons [E > 100 GeV) detected by Fermi with the cat- 
alog of sources detected at lower energies £ < 100 GeV proves 
to be an efficient way to find new y-ray sources observable in the 
VHE y-ray band. Such a study is complementary to the pointed 
observations of selected y-ray sources with the ground-based y- 
ray telescopes. It has enabled production in an all sky survey of 
high Galactic latitude VHE y-ray sky. 

The high Galactic latitude sky survey above 100 GeV has re- 
vealed a large number of new VHE y-ray sources. Addition of 
the new sources to all the sources previously detected with the 
ground based y-ray telescopes doubles the number of known ex- 



tragalactic VHE y-ray sources. Similarly to already known VHE 
y-ray sources, most of the new sources reported in Tables 1-3 
are BL Lac type objects, a special type of AGN with jets closely 
aligned with the line of sight. 

Majority of the sources listed in Table [T] are real new VHE 
y-ray sources, only 2 or 3 sources are expected to be false de- 
tections due to the chance coincidence of arrival directions of 
y-rays with the source position. Based on the analysis of the 
spectral properties of the source in the 0.3-100 GeV band, the 
most probable false detections in Table [T] are source # 33, 34 
and 38. On the other side, several sources in Table 1 have more 
than one photon associated with them. These sources are de- 
tected with highest significance above 100 GeV. Apart from the 
known VHE y-ray sources, the unidentified Fermi source IFGL 
J0505. 9+6121 is detected with significance close to 5cr signif- 
icance above 100 GeV and the sources # 1, 8, 9 and 18 are 
detected with ~ 4cr significance. In addition, BL Lacs with an 
associated front photon whi ch are found to correlate with the 
IVeron-Cettv & VeronI (l2OI0l) catalog (marked by "Vf in Table 
[U should also be considered as high-confidence detections. 

Table[2]lists 21 lower significance sources, among which 5 or 
6 are most probably false detections. Similarly to Table 1, source 
marked "Vf (# 50) which contributes with a front photon to the 
correlation with the BL Lac catalog should be considered as a 
higher significance detection. For the sources listed in this table, 
it is reasonable to accumulate larger exposure to single out false 
detections. 

Table 3 lists VHE y-ray emitting BL Lacs which are not 
listed in the 1-st year Fermi catalog, but contribute to the cor- 
relation of the Fermi VHE y-rays with the BL Lac catalog of 
IVeron-Cettv & VeronI ( |20IO|) . The absence in the Fermi catalog 
might be explained by the very hard spectra of these sources. We 
have demonstrated that 2 of the 7 sources listed in Table 3, the 
sources # 3 and 6 are, in fact, detected with significance ~ 9cr in 
1 -300 GeV band. The source RX J09 1 30-2 1 03, which has 2 pho- 
tons associated with it is individually detected with significance 
more than 4cr above 100 GeV. The other 3 sources which have a 
front photon associated with them should also be considered as 
most probable true detections. 

Some of the new VHE y-ray sources listed in Table [1] is sit- 
uated at large redshifts z ~ 1. Our analysis of the y-ray spec- 
tra of the high redshift sources shows that at least several of 
them do not show any signatures of high-energy cut-off" up to 
300 GeV. Typical fluxes of the sources listed in Tables 1-3 are 
at the level EFe{E > 100 GeV) > (few) x 10"'^ erg/cm^s, i.e 
at the level of 0.1 Crab units. Such sources should be readily 
detectable with the ground based Cherenkov telescopes at the 
energies several hundreds of GeV. Larger collection area of the 
ground-based telescopes should result in the much higher sig- 
nal statistics, which should allow a study of the detailed spectral 
properties of these new bright VHE y-ray sources. We anticipate 
that detailed high sensitivity study of the high-redshift VHE y- 
ray sources with the ground based y-ray telescopes will open a 
possibility to study cosmological evolution of the BL Lac prop- 
erties, of the EB L and of the associated evo l ution of galaxies pro- 
ducing the EB L (iGould & Schrede:^. '1967'; 'Kneiske et all 12004 ; 



Mazin & Raue, 2007; Franc eschini et al 



Stecker etall 120061: 
20 Q8D as well as of cosmological intergalactic mag netic fields 
^Neronov & SemikozLl2009tlNeronov & VovkLl20Ioi) . 
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